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Abstract: We look at the recently proposed idea that susy breaking can be accomplished 
in a meta-stable vacuum. In the context of one of the simplest models (the Seiberg-dual of 
super-QCD), we address the following question: if we look at this theory as it cools from high 
temperature, is it at all possible that we can end up in a susy-breaking meta-stable vacuum? 
To get an idea about the answer, we look at the free energy of the system at high temperature. 
We conclude that the phase-structure of the free-energy as the temperature drops, is indeed 
such that there is a second order phase transition in the direction of the non-susy vacuum at 
a finite T = Tc . On the other hand, the potential barrier in the direction of the susy vacuum 
is there all the way till T ~ 0. 
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1. Introduction 



It has recently been proposed ||l| that susy-breaking in a meta-stable vacuum is an interesting 
paradigm for model-building. The idea is that in field space, far away from the supersym- 
metric vacuum, we could have local minima of the effective potential with non-zero energy. 
These vacua can have parametrically long lifetimes because they are protected by tunneling, 
so using them as phenomenologically viable candidates for susy-breaking is a natural possi- 
bility. Even simple models seem to exhibit the existence of such meta-stable vacua, so it is 
reasonable to hope that this is a rather generic feature of supersymmetric theories. Since 
the original paper came out, many authors have tried to extend this idea in many directions 

1, 0, 1, 1, 0, g, 0, 0, 0, 0, |i8|, ig. 

hi this paper, we will be interested in this question from a cosmological perspective. 
Just the fact that there exists a meta-stable vacuum in the zero-temperature field theory, 
is not enough to guarantee that we will end up in it (and not in the susy- vacuum) , as the 
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Universe cools from high-temperatures. One way to get a better understanding of this issue 
is to look at the phase structure of the free energy at different temperatures. This is precisely 
the purpose of this paper. We work in the context of one of the simplest models presented 
in namely the S-dual of Super QCD^, and look at the theory at finite temperature. The 
scalars in the theory are dual squarks and mesons. At zero-temperature, this theory has a 
supersymmetric vacuum for large values of the meson field. This is the usual susy-vacuum 
of SQCD, as seen from the dual picture. But the interesting thing about the dual picture 
is that here, we can identify new meta-stable susy-breaking vacua of the effective potential 
which are invisible from the pure SQCD picture. This was the key insight of Q. A schematic 



picture of the effective potential at zero temperature is provided in figure 1(c) 



How does this scenario change when we turn on temperature? We can get some idea 
about the situation by looking at the equilibrium thermal properties of this theory at finite 
temperature ||2|, ^, ^, |5|. We could calculate the free energy as a function of the quark and the 
meson fields. At high enough temperature, we expect that free energy has a trough at the 
origin of field space, because the fields are massless there and the entropy is therefore higher. 
We could calculate the mass-matrices for the fields in the quark and meson directions to get 
an idea about the free-energy in those directions. Indeed, when we do this, we end up finding 
that as the temperature is lowered, the second order phase transition in the quark direction 
happens first. There is no (second order) phase transition in the meson direction, all the way 
down to T ~ 0. From these, we conclude that as the temperature drops, the Universe winds 
up in a susy-breaking phase. 

This argument implies that the meta-stable vacuum is plausible, but we should add that 
this in itself is not completely conclusive. The free energy is a purely equilibrium quantity and 
the dynamics of the fields as they interact with a thermal bath could be more complicated. To 
fully understand the situation, we need to do a calculation that incorporates finite temperature 
effects and dynamics. We need to use the evolution of the field at finite temperature: we could 
use the imaginary-time formulation and calculate the friction-type term for the field-equation 
in the thermal bath. We do not report on this calculation here, but the punch-line of the 
preliminary calculation is that at least for some initial values of the fields, we do evolve and 
end up in the susy-breaking vacuum. The free energy equations give us an idea about this 
phase structure, as is clear from the pictures below. In each picture two separate situations 
are plotted together: one in the scalar mesons direction with zero squark vev and another in 
the squark direction with zero scalar meson vev. In l(a)| , at high temperature the free energy 



drives the scalar expectation values to zero. In 1(b) , the plot is at lower temperatures close 



to a threshold where the free energy in the squark direction develops a minimum away from 



zero, while in the meson direction it develops a potential barrier. Finally in 1(c) , at zero 
temperature, the meta-stable vacuum in the squark direction remains cosmologically stable 
O due to the large potential barrier in the scalar meson direction, where the susy vacuum is 



with appropriately chosen number of colors and flavors. 
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located. 



One important assumption is that the Seiberg transition is a post-inflation phenomenon. 
The reheating temperature is higher than the Seiberg transition; from the phase where the 
weakly coupled degrees of freedom are that of SQCD (electric phase) to the phase where the 
weakly coupled degrees of freedom are that of its Seiberg dual (magnetic phase). Besides, the 
unsolved problem of initial value of the scalar mesons after the Seiberg transition allow other 
possibilities that are not considered in the present work. 

If initially, the scalar mesons are located far away in field space, they will execute large 
amplitude oscillations. Then, if both the Hubble friction and the interactions of the meson 
fields with the heat bath are small enough, these oscillations might bring and keep the meson 
fields in the vicinity of the susy vacuum. 

With these assumptions the scalar fields are not in equilibrium with the rest of the 
universe and therefore we should not rely on an equilibrium calculation of the free energy to 
describe this system. It is then plausible, that as the universe cools, the meson fields find 




Q Q 



(a) Effective potential for T (b) Effective potential for T ~ Tc 




Q 



(c) Effective potential for T = 

Figure 1: Evolution of the effective potential witli temperature 
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themselves trapped in the susy state. However in an adiabatic evolution of the universe it is 
appropriate to use thermodynamics of equilibrium to describe the system. This is what the 
body of this paper will make explicit. 

In the next section, we introduce the model under consideration as it was described in 
In section 3, we begin the calculation of the equilibrium free energy of this theory along the 
quark and meson directions by calculating the mass matrices for the appropriate fields and 
expressing the free energy in the high-temperature limit. The interpretation of this object as 
the temperature is lowered is the subject of section 5. We draw some conclusions based on 
our results and end with some possibilities for future work in section 6. 

2. Vacuum Structure of the Magnetic Dual of SQCD 

The meta-stable susy breaking vacuum is in a region of field space which is easier seen in the 
IR (magnetic) dual of SQCD. We report some of the results regarding this theory (following 
the original paper) in this section, but this is by no means an exhaustive review. 

The superpotential for the theory (including the non-perturbative piece arising from 
gaugino condensation) is, 

W = hTi: q^q- hfj,'^ Tr$ + AN{det<^>y/^ (2.1) 

where 

.4 = /l''A^-^+^ v = Nf/N. (2.2) 

The matter content can be described by (the columns denote the SU{N) gauge and global 
symmetry groups): 





SU{N) 


SU{Nf) 


SU{Nf) 


U{1)b 


^(1)' U{l)n 




1 


□ 


□ 





-2 2 


Q 


□ 


□ 


1 


1 


1 


Q 


□ 


1 


□ 


-1 


1 



The second piece in the superpotential breaks the SU{Nf) x SU (Nf) down to a diagonal 
SU{Nf) subgroup. In our notation, the q stand for the (dual) quarks, the $ are gauge singlet 
mesons and is the dynamically generated scale (the scale of the Landau pole) for the IR 
theory. The number of colors of the magnetic theory A^, and the number of flavors A'^^, satisfy 
A^^ > 3N so that the theory is IR free. Notice that all our notation is defined in terms of the 
dual theory and not directly in terms of the microscopic theory (which is SQCD). 



-4- 



For small meson fields, the non-perturbative term can be dropped, and we can calculate 
the moduli space of susy-breaking tree-level vacua. These are at 





^= n ' ' Q= QO), (2.3) 



with QQ = yU^I^r where Q and Q are N x N matrices and (/? is a {Nj — N) x (Nj — N) matrix. 
In this vacum the scalar potential has the value 

Vmin = {Nf-N)\h^^i^\ 

The point of maximum global symmetry in this moduli space, up to gauge transformation 
and flavor rotation, is 

(^ = 0, g = Q = /xIjv, (2.4) 

The interesting point is that along with these meta-stable vacua, far away in the moduli-space 
at large meson values, we have super symmetric vacua that arise from the extremization of 
the superpotential. These are of course the familiar SUSY vacua of SQCD, seen from the 
dual picture. 

Putting all these together, the vacuum structure of the IR dual of SQCD looks schemat- 



ically like Fig. 1(c) . The susy-breaking vacuum is in a local trough (in all directions, even 
though we indicate only the quark direction), so it is meta-stable and it is cosmologically 
stable due to a large potential barrier. 



3. Free Energy 

Our aim now, is to look at the above theory, and calculate the finite temperature effective 
potential (free energy) up to one loop. We want to see how the picture above changes when we 
turn on temperature. This will give us a clue about the possible phase transitions that could 
happen, as the universe cooled. The standard procedure for calculating the finite-temperature 
effective potential is to shift the relevant background fields (in our case the scalar quarks and 
mesons) and use the resulting quadratic pieces in the action to do the computation. So 
essentially, we need to know what the masses^ of the various fields are, as a result of the 
shifts in the backgrounds. It is the calculation of these mass matrices, that we undertake 
next. For the sake of simplicity, we will take the parameters h and /j, as well as the shifts in 
the background fields to be real numbers. 

■^To be more precise, we should say that we are after the coefficients of the pieces in the action that are 
quadratic in fields, we will call them masses even though we are not necessarily looking at the extremum of 
the potential. 
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# of weyl fermions 


fermion fields 


mass 


2N X Nf 


i^q, Ipq 


hif 


1 




A^p^-^iu-l) 


(iV| - 1) 







Table 1: fermions masses, meson direction 



3.1 Mass Matrices in the Meson Directions 

We first calculate the masses of the various fields when a background field is turned on in the 

meson direction, with the scalar quarks set to zero. In the next subsection we will turn on the 
squarks and turn off the mesons. Once we know the mass matrices in both these directions, 
we can draw some conclusions about the free energy and the phase structure as a function of 
temperature. 

Decomposing the meson chiral superfield into a background, trace and a traceless part, 

we can expand the nonperturbative term in the superpotential: 

A7V(det$)V^ = AN^^ex.^ ("^TrMl + + |; 

AiV^'^exp ( _ ^ _ + J 



^From this, we can read off the terms that contribute to the fermion masses in the Lagrangian: 

CD hip T>(V',V9) + -^^^ -i^4>oi^4>o - I^IV'^^^) + h.c. (3.2) 

showing that there are N x Nf Dirac fermions with mass h(p; 1 Major ana fermion with mass 
Aip'^'''^[i' — 1) and {Nj — 1) Majorana fermions with mass Aip^~'^. 

To calculate the scalar masses from the scalar potential, we first compute the F-terms. 



/Nf V y/^f 



hqq + A^'^ ( -4 - KA'PoA (3.4) 



F, = h{4>qY D h{^qf (3.5) 
= h{ct>qf D hiipqf (3.6) 
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^ of real scalars 


scalar fields 


mass"^ 


2{Nf - 1) 


4> 




1 


Im (j)Q 


|yl|2^2.-4(^_l)2_ 

-Re [{A^"-^ - hn'^y{v - - 2)A^''''^] 


1 


Re (po 


Oil A f -.\0 

|^|2(^2i^-4(jy _ 1)2_^ 


N X Nf 


Re{q + q^)/V2 


/iV^ + ReiA^"-^ - h^i'^)h 


N X Nf 


/m(g + g*)/\/2 




N X Nf 


Re{q - q')/V2 


h'^f'^ - ReiAif''-^ - hi?)h 


N X Nf 


Im{q + g*)/V2 





Table 2: real scalar squared masses, meson direction 
Consequently, the scalar potential will have the following quadratic terms: 



Vscalar ^ {Aif' 



u-1 



Tr{qq)h + 



+ C.C.+ 



where the scalar masses are extracted and shown in table 0. 



(3.7) 
(3.8) 



3.2 Mass Matrices in the Quark Directions 

The mass matrices in the quark directions are more complicated than in the meson directions 
because there are contributions from the D-terms. We start by classifying the various sectors 
according to their transformation properties under the global symmetries. We take the shifts 
in the form: (qi) = (qi) = QTL, Q real. The column vectors are N x N^ where Ng = Nf — N. 
The subscript stands for electric. A'e is the same as Nc, the number of colors in the original 
(microscopic) theory, namely SQCD, but we prefer to think of it purely in terms of the dual 
theory. On a similar note, is actually Nm, with the subscript standing for magnetic. 

There are four sectors that do not mix with each other in the Lagrangian. We will use 
this fact to our advantage in calculating the mass matrices. These sectors are: 

1. NeXNe-. (j)22, 

2. NeX Nm ■■ q2,4>21, 

3. NmX Ne : q2,4>12, 

4. NmX Nm ■■ qi,qi,(l)u,V. 

The V in the last line is the vector superfield, it gets massive through a Higgs mechanism. 
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Now we look at the masses of the various fields sector by sector. To start off, in the 
A'^e X A^e sector, there are N"^ x (complex scalars + Weyl Fermions), and all of them remain 
entirely massless (at tree level). But it is in this sector that supersymmetry is broken by a 
positive contribution in the scalar potential. Decomposing ^22 into trace and traceless parts, 
we get 

$0 



$22 = ^2 + $22; F^o^ D -hiiyjNe. 

Hence 



Vscalar > N^h^ fl"^ . 

It is easiest to deal with the two mixed electric/magnetic sectors (2 and 3) together. 
Together, there are 2NeNm x (2 complex scalars + 1 Dirac fermion). The fermionic masses 
arise from terms like 

hQ Tri^f^^l) + hQ Tr(^'f ^^^) + h.c, (3.9) 

and these give rise to a Dirac mass of hQ. To calculate the bosonic masses, we need the scalar 
potential in the quark direction which can be calculated easily enough from the superpotential 
and the F-terms. The F-terms Fq-^ , Fq-^ do not give rise to scalar masses in sectors 2 and 3 
because the scalars from these sectors have zero vevs. The relevant non-vanishing ones are 

F^,, = h{q2qiY 3 /iQ(g2)* (3.10) 

F^,, = Kqiq2)' D hQ{q2)' (3.11) 

i^</.22 = h{Q2q2f -h^?xl (3.12) 

Fq2 = Hqi4'12 + g'2022) ^ hQcpu (3.13) 

Fq^ = h{(j)2iqi + </'22g2) 3 hQ(p2i (3.14) 
Therefore, the scalar potential contains the quadratic terms: 

Vscalar ^ {Q^m^ + Q\2\^ " Tv{q2q2 + h.C.) + Q^l^ul^ + Q^\cP2l\^) (3.15) 

where in our notation modulus squared of matrices means trace over the product of the matrix 
and its adjoint. As we see, 2Nm x A'e complex scalars (</>i2 and 4>2i) get squared mass h?Q'^, 
2Nm X A'e real scalars {Re{q2 ± ql^/y/2) split their masses into h?{Q'^ =b /i^) and another 
2Nm X A^'e real scalars {Im{q2 ± (^^I\f2) get mass /i^Q^. See table (P). 

Now we turn to sector 4. First, we separate out the background: 

gi = QI + ^i ; gi = QI + ii (3.16) 
Some of the fermion masses arise from the terms 

g^Q Tr(A^?) - g^Q Tr(A^'f ) + h.c. (3.17) 
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^ of real scalars 


scalar fields 


mass^ 


^ ^ ' 771 ' 6 


(t>A9 and 1^91 

Y^iZ "'^J-'-J- y^ZL 






Re{q2 + g1)/^/2 




Nm X iVe 


iie(<?2 - g1)/\/2 




A^m X Ne 


/m(^2 + ql)/V2 




Nm X iVe 


Im{q2 - ct2)/V2 





Table 3: real scalar squared masses, sectors 2 and 3 



# fields 


fields 


mass 


iV2 - 1 




2gQ 


Af2 _ 1 
m 


A" 


2gQ 


- 1 


traceless - ^?)/\/2 


2gQ 


2N'^ 




hQV2 




(^? + M/?)/V2 


hQV2 



Table 4: vector boson and fermion masses, sector 4 

where gauginos A*^, and the traceless (^'^ — ^'^)/\/2 have equal masses 2gQ. 

^From the Kahler potential, the vector bosons get a mass 2gQ. The traceless part of 
Im (^{qi — qi)/\/2^ is gauged away and the traceless part of the scalars Re (^{qi — qi)/\/2j 
get their masses from the F-terms (to be calculated) and from the D-terms. The contribution 
to this squared mass from the D-terms can be read off from 



9 



.^(QTr(2i?er(gi = ^ (Tr(i?er (gi - |i)/V2 



(3.18) 



to be Ag'^Q'^. The scalar /m(Tr(gi — qi)) /V2 and the fermion Tr(^'^ — ^'^)/\/2 remain 
massless at tree level. But i?e(Tr(gi — qi))/\/^ receives a contribution from the F-terms. 

The shifts in qi and qi give rise to terms of the form 

hQ^f^i^l + ^l) + h.c. 

From these, the 2N^ Weyl fermions (^-f^) and (^'f + ^'f)/\/2 acquire a mass of hQ^/2 each. 
Their respective scalar superpartners acquire mass through the scalar potential. Writing the 
relevant terms in the F-terms, using ( p. 161 ), we have 
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# of real bosons 


fields 






Re {qi +'ii)/V2 


/l2(3Q2_^2) 








Ar2 

^ ' m 


Im{qi +'g~i)/V2 




_ 1 


traceless Re{qi — qi)/\/2 


/l2(;,2_g2)+4^2Q2 


1 


Tr(i?e(gi-|i)/y2) 


/i2(;x2-Q2) 


1 


TV(lm(gi-|i)/V2) 






Table 5: real squared masses for bosons, sector 4 



and D hQ(pii , i^'g^ D hQ(pii. Hence, the quadratic terms in the scalar potential are. 



Vscalar ^ h"" I 2Q2 | 1 2 + 



Qi + qi 



V2 



+ (Q2_^2W Tr(gi gi) + c.c 



(3.19) 



This shows that among the real scalars, 2N^ get squared masses 2/i^Q^, N!^ get squared 
masses h?{'iQ'^ — /i^), and N"^ get squared masses 2/i^Q^. The corresponding fields are ^n, 
Re{qi + qi)/y/2) and Im{qi + qi)/\/2 respectively. The terms Re{qi — qi)/^/2 get mass from 
above and from the D-terms, splitting the field matrix into a trace part with mass h?{^'^ — Q"^) 
and N'^ — 1 traceless components with mass — Q"^) + Ag^Q'^. 

The vector boson, gaugino and fermion masses are presented in table (Q), and the scalar 
masses are in table (0). 



3.3 Effective potential 

The effective potential at finite temperature is the free energy of a system in a thermal bath. 
In thermal equilibrium, there is no time dependence and the different phases correspond to 
local minima of the free energy. 

The finite temperature effective potential up to one loop is given by Q 



where Vtreei'Pci) is the classical piece. The one-loop correction, {(j)ci), for a generic theory 
is: 



Nb + -Np H 

90 \ 8 / 24 



T 

12^ 



j;(M|), + 3j;(M2), + j;(M|) 

i 



+ 



(3.20) 



- 10 - 



and Vi{(j)ci) is the zero-temperature piece, {Ms)i, (My)j and {Mp)i rae mass-matrix eigen- 
values for the real scalars, vectors and Weyl fermions respectively, and A'^ = Np is the 
number of bosonic/fermionic degrees of freedom, paired by supersymmetry. What we have 
done here is to follow the standard practice and split off the one-loop, finite-temperature 
effective potential into a part that is independent of temperature (and therefore is the same 
as the zero-temperature effective potential) and then do a high-temperature (T 3> masses) 
expansion on the remaining (temperature-dependant) piece. The zero temperature piece is 
calculated using the usual supersymmetric generalization of the Coleman- Weinberg formula 
§: 

1 AA"^ 

where A4 stands for the full mass-matrix, with the supertrace providing the negative sign for 
the fermionic terms. 

Armed with the above expressions and the mass-matrices from the last section, its easy to 
calculate the finite-temperature effective potential. In the following, we only keep the terms 
that are quartic and quadratic in temperature. In the case of background fields in the meson 
direction, we find: 

Vtreei^) = iV/K^C^""' " hfl'')\\ (3.22) 

.0,,, 1 (_,|.|V-(.-i)'.„gl^lV-V-i)V 



64^2 V -i^^i ^ A2 
+ [\A\^ip'^''-^{u - if - Re [{Aip''-^ - hi?y{v - - 2)Aip''-^]fx 

Ml V^''"^(i^ - 1)^ - Re [(^93^"^ - h^i'^Yiv - l)(v - 2)A^^-^ 
X log ^ p 

+ [\A\'^^'^''-^{v - if + Re [{A^''-^ - hi?y{v - l){v - 2)A^''-^]fx 
\A\'^if'^''-\v - If + Re [{Aif'^-^ - h^'^Yiv - l){v - 2)Aip''~^] 



X log 



A2 



+N^Nfh\h^^ + Re{A^^-^ - h,^)f log h^^ + MA^I^^jH^ _ 

-2N^NjhV log ^ + N^Nfh\h^' - Re{A^^-^ - h^')f log - " ^^')^ 

(3.23) 

Vn^) = -^{{Nf + N^f-l) + 
+^[3(iV| - 1)|A| + 3|A|V''"^('^ - 1)' + 8N^Nfh^ip\ (3.24) 
Similarly, for the quark direction: 

VtreeiQ) = N^h^ fi^ + N^h\Q^ - fi^f, (3.25) 
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ViiQ) = ^{-2NM^Q' log ^ + NM\Q' + f^'? log + 



+Ar^iVe/^^(g2 - log + iV^/.4(3Q^ - „r log " 



^2 ' m'" V"^ > ^2 

-4iV2 /,4q4 log ^ + (iV^ - l)(/,2(^2 _ q2) + 4^2q2^2 W^^+W ^ 

- Q^)f log _ (iV^ _ 1)(VQ2)2 log 1^), (3.26) 

yI{Ql) = ((^/ + ^^f - 1) + ^ [3iV^^e/i' + l^lh" + 5' + 5(iV^ - 1)5'] . 

(3.27) 

With these explicit forms for the finite temperature effective potential, we will be able to 
draw some conclusions about the nature of the phase-transitions in the next section. 

Armed with these expressions and the mass-matrices from the previous seb-sections, we 
can calculate the explicit forms for the free energy. 



4. Cooling and the Emergence of Different Phases 



We want now to understand what happens to the mesons and squarks 9?, Q during the evo- 
lution of the universe, as we cool down from a high temperature. In particular, we want to 
know the phase structure. If the phase transition in the quark direction happens at a higher 
temperature than in the meson direction, then we have at least some reason to believe that 
we will eventually end up in the susy-breaking phase. 

Phase transitions are characterized by a critical temperature Tg. By definition, the critical 
temperature Tc for a second order phase transition is the temperature at which the second 
derivative of V^(</?, T) at the origin, in one of the field directions, changes sign from positive 
to negative. When this happens, the local minimum (at the origin) in that direction becomes 
a local maximum and a new minimum forms at some finite field value. As a consequence, the 
vacuum at the origin becomes unstable, a phase transition takes place, and the fields evolve to 
the newly formed minimum. Of course, again, we emphasize that we are doing an equilibrium 
analysis, but we believe that this is enough to give a preliminary, heuristic picture of the field 
history. We find that in the quark direction is given by 



On the other hand, in the meson direction, we see that there is a minimum away from 
the origin at this temperature, but, the local minimum at the origin is still there. At some 
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temperature Tc these two minima will become degenerate {V{0, Tc) = V{ipm-,Tc)), but there 
still is a potential barrier between them. Tunneling through the barrier can start when the 
temperature hits Tc , but this phase transition is first order as opposed to the second-order 
phase transition in the quark direction. The critical temperature for the first order phase 
transition in the mesons direction turns out to be 



We notice that in the meson direction the origin is always a local minimum for every 
temperature, even zero temperature. In fact expanding the 1-loop effective potential at zero 
temperature V{ip, Q,T = 0) around the origin we find 

V{^,Q,T = 0) ~ /lV^ ~ 0. 

Finite temperature effects do not change the fact the the origin is still a local minimum in 
the meson direction. In this case we have 

V{ip,Q,T) ~ T^h^ip'^, ~ 0. 

Because the phase transition in the meson directions is first order, it is accomplished through 
quantum tunneling processes and hence is much more strongly suppressed than the classical 
phase transition in the quark directions. 

To gain a better understanding of the phase transition in the quarks direction we studied 
the finite temperature effective potential for every value of (j) and Q close to the critical 
temperature . The result of this analysis is plotted in Fig. |2| From the shape of the 
effective potential around the origin we immediately realize that the flow of the vev happen 
in the Q direction and it is not possible for the vev to flow in the (j) direction. 

We are now in the position to form an idea about the phase history as the universe 
cools down. Let's suppose that we are starting at a temperature T 3> Tc. We could for 
example be in the reheating phase after inflation. At this temperature, the origin of field 



space is a minimum for the finite temperature effective potential V{if, Q, T), figure 1(a). This 
is qualitatively plausible, since the massless fields make the biggest contribution to entropy. 
We also make the assumption that when we start off at this high temperature, the mesons ip 
and the quarks Q are localized around the origin of field space: 93 = 0, Q = when T ^ Tc 
As the temperature decreases to T = T^ , the curvature of the effective potential V{(p, Q, T) 
at the origin becomes negative in the Q direction but it stays positive in the direction: 



dQ 



<0, 

<^=0,Q=0,T=Tc 
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1.5 ° 

Figure 2: Effective potential for every values of and Q for T ^ 

'd^v\ 

> 0. 

</p=o,(3=o,T=rc 



dip 



Also, at T = T^, a new minimum forms in the Q direction, see figure 1(b) . As 



a 



consequence, a phase transition occurs and the fields move to the newly formed minimum 
Qm- As the temperature of the universe continues to decrease, we eventually arrive at T ~ 0, 



see figure 1(c), and the minimum Qm becomes the (meta-stable) non-supersymmetric vacuum 



= j ■ meson side, as the temperature drops, the minimum ipm becomes 

the supersymmetric vacuum ip^ = ^ (Nf-3N)/{Nf-N) i but thankfully, phase-transition into 
the susy phase is suppressed by tunneling at all stages. In writing the expression for the 
susy- vacuum ip^, we use the Intriligator et al. convention, with e = /x/A^ where Am is the 
dynamically generated scale of our (infrared) theory. As mentioned earlier, it is the scale of 
the Landau pole. 

Thus, the phase structure of the theory seems to imply that for reasonably tame initial 
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conditions for the scalar quarks and mesons (namely, they start off near the origin of field 
space), the phase transitions lead us into the susy-breaking vacuum at T ~ 0: 



5. Conclusions and Future Directions 

We found that the phase-structure of the free energy seems to imply that as the universe cools, 
we end up in the susy-breaking vacuum ^ . In this final section, we point out some caveats and 
limitations inherent in our study. One of these, is the assumption of thermal equilibrium. By 
working with the free-energy, we are ignoring the possibility that the fields could evolve and 
interact with the heat bath. The dynamics of the fields could result in overshoot or undershoot 
around the vacua. But it seems unlikely that these effects will destroy our conclusion if the 
scalars are starting at the origin of field space, at least if the evolution of the universe is close 
to adiabatic. That brings us to a related issue, which is that we don't have a good idea about 
what are good initial conditions for the scalars, as was discussed in the introduction. 

When calculating the effective potential, we focussed on just the meson and squark di- 
rections. There is the possibility that we would discover new valleys and slopes if we were to 
make a map of the potential for all field values. One possibility is to do a perturbative analysis 
around the quark/meson directions to see whether these directions are indeed troughs and 
not saddles at all temperatures. 

Another issue which could affect our conclusions, even though it is likely to be less 
important, is the question of the time scales of first order versus second order transitions. It 
seems likely that the original argument of Intriligator et al. suggesting that the meta-stable 
vacuum can be made parametrically stable against tunneling, should work (perhaps with a 
suitable modification) for our case too. 

To deal with non-equilibrium situations and initial conditions, one could work in the 
context of the real-time formulation of finite temperature field theory where we can consider 
both dynamics and thermal effects. There, the interaction with the heat bath would be 
encoded in a friction-like term. Preliminary investigations in this direction seems to imply 
that at least for some choice of scalar initial conditions, we end up in the meta-stable vacuum. 
We leave a more thorough analysis of this issue for future work. 
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